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Введение. Исследование проникновения упругих волн через 
периодические решетки является важной проблемой в обла- 
сти ультразвуковой количественной оценки материалов, рас- 
пространения звука и для электромагнитных волноводов с 
диафрагмами. На практике аналитические результаты могут 
быть получены в предположении низкой частоты, со слабым 
режимом взаимодействия, когда лишь некоторые прибли- 
женные результаты можно установить в аналитической фор- 
ме. 

Материалы и методы. В предыдущих работах автором 
изучены 3-О задача проникновения волны нормальной 
плоскости с бесконечной двоякопериодической системой 
трещин в низкочастотном режиме и 2-0) задача 
проникновения волны нормально двум системам, когда 
каждая содержит бесконечный периодический массив 
трещин.Целью настоящей работы является обобщение 
полученных ранее данных —Й результатов исследования 
свойств рассматриваемой системы, основанного на плоской 
задаче о распространении волн в упругих средах с 
периодическим массивом трещин. 

Результаты исследования. Настоящая работа продолжает 
изучение 2-) задачи для трех параллельных массивов, обра- 
зующих двоякопериодическую систему. Исследование по- 
священо выводу аналитических выражений коэффициентов 
отражения и прохождения, когда плоская продольная волна 
падает на систему трех идентичных плоских решеток, распо- 
ложенных друг за другом. В режиме частотного диапазона 
одной моды задача сводится к системе гиперсингулярных 
интегральных уравнений, решение которой дает эти коэффи- 
циенты и явное представление волнового поля внутри струк- 
туры. 

Обсуждения и заключения. Применяемый метод позволяет 
управлять акустическим фильтром в рассматриваемой ре- 
шетке выбором соответствующей длины трещины, частоты и 
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Введение. Исследование проникновения упругих волн через периодические решетки является важной проблемой 
в области ультразвуковой количественной оценки материалов, распространении звука и для электромагнитных волново- 
дов с диафрагмами. Различные численные методы были применены в двумерных задачах с периодическими отверстиями 
произвольной формы [1-3]. Результаты компьютерного моделирования продемонстрировали высокую точность. При этом 
развиваются и аналитические теории. На практике некоторые приближенные результаты можно установить в аналитиче- 
ской форме в предположении режима низких частот при слабом взаимодействии волн. Таким образом, аналитические 
методы, обеспечивающие явные формулы для соответствующих параметров рассеяния, задают только определенный низ- 
кочастотный предел. В работах [4—7] получены явные аналитические формулы: 
— для параметров отражения и прохождения в режиме одной моды акустической волны, проникающей сквозь двояко-, 
троякопериодические массивы отверстий и объемных препятствий произвольной формы; 
— для плоских задач распространения волн через периодический массив трещин в упругих твердых телах [8]. 

Случаи высокочастотной дифракции на трещине анализировались в [9-11]. 

В данной работе изучаются три массива, каждый из которых представляет собой периодическую систему трещин 
в двумерной постановке. Как показано в [12], качественные свойства такой геометрии аналогичны соответствующей 
плоской задаче для двух массивов. Предполагается, что при нормальном падении волны имеет место режим распростра- 


о © Ш х № 
нения одной моды, т.е. а. <л, где К, — волновое число поперечной волны е ", 2а — период данной системы. А 


также расстояние 4 является таким, что выполняется условие 4 /а>>1. 

Целью настоящей работы является обобщение результатов, полученных ранее в [12, 13], и вывод новых аналити- 
ческих выражений коэффициентов для системы трех параллельных экранов (с учетом свойств ядра гиперсингулярного 
интегрального уравнения). Соответствующие вычисления основываются на решении плоской задачи распространения 
волны в упругом теле с периодическим массивом трещин. 


Постановка задачи. — Вывод интегральных уравнений. Рассмотрим структуру, состоящую из трех бесконеч- 
ных массивов коллинеарных трещин, вертикально расположенных в х = 0; 4; 24. Каждая из трещин, длина которой 2, 


о < Хх 
расположена с периодичностью 2а по оси у. При изучении падения плоской продольной волны е Г на данную систему 
вдоль положительного направления оси Хх, в силу симметрии вопрос СВОДИТСЯ К рассмотрению волновода ширины 2а 
(рис. 1). 




















Рис. 1. Фрагмент периодического массива трещин 


Ес. 1. Зестеш оРрепо(1с аггау оЁ сгаск$ 
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Тогда потенциалы, удовлетворяющие уравнению Гельмгольца, есть функции: 


п=1 


И Е № | 
фи" = ем" + Ве "+ У Але" со 
а 
щие = _ ет" ‘(2 х<0; 
а 


фе" + Е соз [А (хх, + НЕ СОЗ (хфх,)]+ У4Е*с Ма, (хх, + НсНа, (хх Зрекк[ у } 


п= 


ли’ = Убе (хх, 1] Рус, (х—х, оз а] (5—Па <х<54; 5=12; 
п=1 а 


;: Ш © —а [х-2а 
фе — Тейа-24] УСе 4х 1вих( 
п=1 а 


це = Уре" ни [ х>24. (1) 
а 


п=1 


Заглавные буквы здесь являются неизвестными константами, а 


4, = [я / а)’ ы = [и / а)’ -в |. (2) 


Гармонический временной множитель берется в виде е '°', К, К, являются продольным и поперечным волновы- 
12/2 


ми числами, а К, Т — коэффициенты отражения и прохождения соответственно. Далее возьмем компоненты тензора 


напряжений о,.; 0, и вектора перемещений и,,и, в терминах потенциалов, используя стандартные формулы Грина — 


ж? Эх 


х Жх 
Ламе. В рассматриваемой задаче плоская продольная волна с потенциалами фу =е "Г, \, =0 приходит из —0 и порож- 
дает отраженное поле. Предполагая непрерывность поля перемещений вне трещин, введем следующий вектор неизвест- 


ных функций &,(у), 8,(у), |УКЬ: 
и и = 21(), х=0; пФ -и® = 27), х=а; пФ -и = Бу), х=2а. (3) 
Используя (1)-(3), найдем все необходимых постоянные, входящие в потенциалы (1), через &,(у), &,(у). Для 


этого сначала проинтегрируем уравнения (3) по отрезк <а, чтобы получить систему для 6 коэффициентов, не 
У у м» 


входящих под знак суммы В выражениях для потенциалов: 
ГВ 
—В- НЮ за ва) = — [2104 ; 
2а -Ь 
1 . 2 : . 2 
—Р А зт(а)- Но зщ(а) = 5а | 504; 
-Б 


Ь 
РАЗ зи а) + ем“ ит = -- [а(4». (4) 
а 


Свойство ортогональности тригонометрических функций с учетом предположений метода 4,4 >>1, п=1,2,3... 


приведет (3) к соотношениям, связывающим остальные 12 о 


[А+ НоэИ(а, У» ГВ, — Рен = 5 р ато] п ш: 
а 


[-4, + лева а) В, + РАЗ), 5} а Е 2 
а 





А+ НУ, 54,4) +6, — Ни) = ры 
2а —ь а 


Ь 
[21+ 52] “са, а) [61 + Р?]и. (4) = =] 820) ни [ 2%) ф; 
а а- а 


Ремизов М.Ю. Распространение упругих волн через периодическую систему трещин 








Ь 
ГЕ25(а„а) + Са, 2.) - Б,]^^ = >. ТВ. рии "6 
а а 


[-Ерсй(4„4) + с," —[С,зй(,4) +2, |, =! 807) эт 7 "|. (5) 
Условия непрерывности напряжений на границе имеют место в виде 
0; 00 =с@, х=4; с® = 00; о = 000, х=24. (6) 


52 = 0; в“ = в, х=0; 00) =с@); г 
После применения для (6) аналогичной процедуры интегрирования по отрезку и принимая во внимание свойство 
ортогональности тригонометрической функции, получим остальные 9 соотношений, которые вместе с (4) и (5) образуют 
алгебраическую систему с 18 неизвестными постоянными задачи. Часть решения системы, соответствующая х=0, 


Н®,Е®,Н ый К ‚ выражения которых, используемые ниже, принимают вид: 
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2 2 
| 1 7 + ап Ш Ч нь М 
п [ С дао 5 а | „(п)зт| —— |4]; 

2а5й(4„4) «8 а ь а 1 а 


2+ 
А Вы ие тов" . пт — 1 азодза[ 71 - пит 

















” 2 аъ 5 -Ь 
12а 
В ет [п пт | вродь Пит 
ра 1 ю -Ь а ср, т -Б 
1 Ц 
В = ых т без ИИ "а ( он 1. (7) 
аа) = и й а р с, т тр ы | 


Формулы для остальных постоянных для краткости опускаем. 
Теперь полученные постоянные используются для удовлетворения условий равенства нулю компонент тензора 


напряжений на трещине. Здесь следует принять во внимание, что в задаче остаются только нетривиальные функции рас- 
крытия 8°(у);| ук; 5 =1,2,3 ‚ учитываемые в соотношениях 


(6) (0) о (2) 
Гери ох =0, х=0; ох, =0%. =0, х=а; 0; 
Опуская промежуточные преобразования, окончательно получаем следующую систему трех интегральных урав- 








2 = 009 =0, х=24. (8) 








нений относительно функций 2*(у); [у |< 5; 5=1,2,3: 




















е4 424 ь 
Ком По ег 85а! =1; 
т 

ыы Ей Г. ыы 

р 5. К(’-Ве?(Ва+ (4 =е"; 
дат им - (7-9 82(04 т 1 
Ж 24 Ь ай И 
т т керио -5_ | ко -деш = ет”, (9) 

2 


где К) +У (В, 14,2) соз(лиу а), В, а а, = пп /а. 
п=1 


Свойства интегрального уравнения. Начнем с рассмотрения вспомогательного интегрального уравнения 
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16 К, 
5 Гит)К (у -п)4ч=1; К(у)= С, соз(а,у); Г, = п [У[5Ь. (10) 
-Ь 


2 
Иа 2 Ч 


Заметим, что Г, = —2(1 < / с? )а„, п —> 0. Тогда выражение ядра преобразуется к виду: 


2 2 


2 \х со 
К(>) = 3-5 Хан соз(а„у) + У Ё ‚з- | ввыя К(>) --3 (У) +К,(у) * 


С] п=1 п=1 С] С] 


Здесь вторая сумма ядра есть некоторая регулярная функция. Первая имеет как регулярную, так и сингулярную 
части: /(у) = [1,.(у)+1,(»)], для выделения которых рассмотрим выражение 


“= У псоз( у). (1) 
п а а 


Используя обобщенное представление для суммы (11), имеем преобразование 

2 
5 ИИД ие И" 1 а 
У псо8(— у) = Ши Уе “"исо(— у) = — ‚ С-5, У—>0), после чего ядро вспомогательного 
п— а &—>0 и=1 а 4зт“(пу/2а) пу 





интегрального уравнения (10) представляется в виде: К(у)=-2И- с / с? УГ.) + Г, (›)]+К,(2), где сингулярная и 


а Ил 








регулярная части есть соответственно /,(у) = ы = ЕЕ . 
пу пу 4азш“(пу/2а) 

Легко видеть, что особенность ядра для малых аргументов содержит гиперсингулярность по одной переменной, 
возникающую в теории трещин, хорошо известную в задачах теории упругости в неограниченных средах [14]. Для 
устойчивости численного метода при выполнении расчетных экспериментов здесь применяется дискретная квадратурная 
формула для гиперсингулярных одномерных ядер, известная как «метод дискретных вихрей» [15]. Используя новое 
представление ядра основного интегрального уравнения (10), перепишем его следующим образом: 











2 
с 
ГВ ЩЕ с 
— [11| Ф,(у-п)+ |4 =1, Ф,()=-21-—,0)+К, О). (12) 
2а -ь п(у-п) с 
Дискретизация интегрального соотношения (12) приводит к формуле 
2 2 
241 -“? 241-° 
| м* С] 1% Ц С] 
=>, ] ВФ. Ол) 5 4п= -— Хй(ить) [| #(1)| Ф, С, — пк, + = 4(т-у,) |= 
2а КН пу —п) 2а к= ПА п(у —1) 
2 2 
24-2]  241-“ 
м с 


т 
= — Уи») Ф.П › 
2ака п-т)  пО- 


где |; =-Б+А&,, у =-Ь+(1-0,5)81 , [=1,...,М, & =2Ь/М. 
В игоге имеем алгебраическую систему уравнений: 
2 
с 
2ай1-—] 


гм с: 1 1 
— УЙ(ць)| Ф.О - кв + =1. 
2а ка д (9-м) О-тМьа) 





Аналитические выражения волновых характеристик. Используя вспомогательное уравнение, входящее в (10), 
для решения которого вводится постоянная 


Н= | паи, (13) 
-Ь 


в терминах функции Й(у) из (9) и (10) получаем: 
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а 24 


5? ву (и) = АСТ /4ай ).Л + (е 


5? в (п) = ПСА / 4ай, Л, + (е 


Иа Л, + (ет 14а Л, -П; 


И дик АЕ 29 |: 


223 (п) = РГА / 4а ).75 + (е® /4аК,).Л› + (2? 4, )Л -е “1, (14) 
где обозначено 
ь 
= Г еоаь, у=12,3. (15) 
-ь 


После интегрирования уравнения (14) по множеству | у|<Б с учетом (13) получаем систему линейных алгебраи- 


ческих уравнений относительно „Л р 1=12,3: 


(2 -Н/4ай,)Л ей /4ай Л, — ый > М: =-Н; 
сне" дай) +6? - Н 14а —(е "4 Н дай ль = Не" ; 
(Нем? уда —(е“ Н дай Л, +6? — Н /4ай Л = ие (16) 


После того как найдено решение системы (16), т. е. найдены величины „Л р 1=1,2,3 и параметр Н ‚ находим все 


необходимые волновые характеристики. В частности, в отражения и преломления будут выражены в виде: 

















р. 2“ 24а Ь 
3 
= ВЕ; 
ее Ти > М (па 
24 йа ь 
__е 1 #24 
= Газа. С т 1 &.(чет”. (17) 

Чай —ь Аа 


Отметим, что данная теория позволяет найти ое (17) в явном виде, а единственным вопросом, кото- 
рый решается здесь численно, является построение решения основного вспомогательного гиперсингулярного интеграль- 
ного уравнения (13). Исходя из этого, предлагаемый метод может быть назван «полуаналитическим». 

Обсуждение и результаты. На рис. 2—5 продемонстрированы некоторые примеры зависимости физических 
параметров от частотного параметра а/п, (К =К,) и длины трещины в упругой среде с отношением упругих 
скоростей с, /с› =1,87. Модуль коэффициента прохождения |7 | с ростом частоты после нескольких осцилляций всегда 


уменьшается, независимо от геометрии решетки. Осциллирующий режим и одновременно определенные потери 
прохождения очевидны при увеличении длины трещины Б/а (рис. 2). 








О 0.2 0.4 0.5 0.8 1 


Рис. 2. Зависимость коэффициента прохождения от частоты: а =1; 4 =4; 
кривая 1 —В/а=0,15; кривая 2 — Б/а=0,2 
Ею. 2. Оереп4епсе оРтапзги1 оп соеййстепе оп Недиепсу: а =: 4 =4; 
сигуе 1 — Б/а= 0,15; сигме2 — Б/а=0,2 


Вестник Донского государственного технического университета 2017, №1(88), 18-27 








При фиксированных ак. /л и Я модуль коэффициента отражения | | при изменении параметра Б/а может 


иметь различный характер поведения (рис. 3). 
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Рис. 3. Зависимость коэффициента отражения от длины трещины: 
а=1; 4 =4; кривая 1 — ак, = 0,222; кривая 2 — ак, = 0,357 ; кривая 3 — ак. /л = 0,595 
Е. 3. Ререпдепсе оЁтеЙесноп соеРйслещ оп сгасК?’$ 1еп2 в: 


а=1; 4 =4; сагуе 1 — ак, = 0,222; сигуе 2 — аК› =0,357 ; сагуе 3 — а. /л = 0,595 


Для одних частот это монотонно возрастающая функция на всем диапазоне изменения величины трещины 
(кривая 2). В других случаях при некотором размере трещины достигается локальный минимум (кривая 1) или даже 
несколько локальных минимумов (как на кривой 3). При некоторых размерах трещины может происходить полное 
запирание канала, когда |Т |=1. Если 4 и ак. /п зафиксированы, то максимальное отражение наступает для наибольшей 


трещины в материале. И наоборот, слабое отражение происходит для небольших трещин, что вполне естественно с 
физической точки зрения. 
Также отметим, что для заданных В/а и ак, /л величина | Т(4)| может достигать минимума. При этом имеются 


такие значения 4, для которых достигается полное прохождение (рис. 4). 
7 
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Рис. 4. Зависимость коэффициента прохождения от расстояния 4; Б/а=0,45; ак, =1,87 


Е. 4. Оереп4епсе о гап$и11з310п соеЁйслетЕ оп 415апсе 4; Б/а=0,45; ак, =187 
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На рис. э показано, что свойство акустических фильтров достигается не только для относительно длинных 
трещин, что вполне естественно, но также и при средней их длине. 














Рис. 5. Зависимость коэффициента отражения от частоты: а =1; 4 = 4; 


Кривая 1 — Р/а=0,7; кривая 2 — Р/а=0,9; кривая 3 — Б/а=1 
Ею. 5. Оерепдепсе оРтеЙесноп соеЁЙс1еп оп Недиепсу: а =1; 4 =4; 


сигуе 1-—В/а=0,7; сигуе 2 — Б/а=0,9; сагуе 3 — Б/а=1 


Кривая 3 для максимально длинной трещины подтверждает правильность настоящего расчета при анализе 
предельного случая, когда значение функции отражения становится постоянным почти на всем одномодовом частотном 
интервале, что вполне естественно. В самом деле, для любой геометрии в случае крайне низких частот коэффициент 
отражения начинает расти от нулевого значения. С другой стороны, для любой фиксированной частоты, | Ё |->1, когда 
Ь/а—>1. Поэтому для больших трещин возрастание | А | от нуля до единицы должно происходить в начале одномодоваго 
интервала по частоте. Однако интересно также, что уже для среднего размера трещины (кривая 1, Б/а=0,7;) имеется 
пара относительно длительных частотных интервалов, где коэффициент отражения приближается к единице, при этом 
существуют такие частоты, при которых запирание волны резко меняется на почти полное ее пропускание. Вообще же 
свойство фильтра на верхних частотах данного одномодового интервала имеет место для всех представленных значений 
длины Б/а. 

Выводы. Специфика внутренней структуры в рассмотренной задаче позволяет достичь свойств акустических 
фильтров за счет выбора частоты, расстояния между периодическими массивами трещин, а также длины имеющейся здесь 
периодической неоднородности. Метод, развитый в настоящей работе, может эффективно применяться в исследовании и 
более сложных трехмерных волновых задач, когда вертикальный массив, содержащий двухпериодическую систему 
трещин, последовательно повторяется периодически конечное число раз вдоль направления распространения волны. 
Троякопериодическая система будет предметом следующей работы, в которой предполагается рассмотреть некоторые 
новые специфические свойства волн в упругих средах, не наблюдавшиеся при исследовании аналогичных задач в 
скалярной постановке. 

Автор выражает благодарность заведующему кафедрой теоретической и компьютерной гидроаэродинамики 
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